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We investigate partially disordered antiferromagnetism in CoCl2-2SC(NH2)2, in which a-b plane 
hexagonal layers are staggered along the c-axis, rather than stacked. A robust 1/3 state forms in 
applied magnetic fields which the spins are locked, varying neither as a function of temperature nor 
held. By contrast, in zero held, partial antiferromagnetic order occurs, in which free spins are avail¬ 
able to create a Curie-like magnetic susceptibility. We report measurements of the crystallographic 
structure, and the specihc heat, magnetization, and electric polarization down to T = 50 mK and 
up to fj,oH = 60 T. The Co^^ S = 3/2 spins are Ising-like and form distorted hexagonal layers. The 
Ising energy scale is well separated from the magnetic exchange, and both energy scales are accessi¬ 
ble to the measurements allowing us to cleanly parameterize them. In transverse helds, a quantum 
Ising phase transition can be observed at 2 T. Finally we hnd that magnetic exchange striction in¬ 
duces changes in the electric polarization up to 3 pC/vc? and single-ion magnetic anisotropy ehects 
induce a much a larger electric polarization change of 300 /rC/m^. 

PACS numbers: Valid PACS appear here 


I. INTRODUCTION 

The puzzle of frustrated magnetic spins on a triangu¬ 
lar lattice has intrigued the scientific community for more 
than half a century. [1] Classical Ising spins can exist in 
only two states, and so requiring them to accommodate 
a three-fold lattice with antiferromagnetic interactions 
is a rich source of interesting magnetic patterns on the 
micro- and meso-scale. The ’up up down’ or ’up up up’ 
arrangements of Ising spins in a triangle provides the 
smallest magnetic unit cells. Achieving a continuously 
varying magnetization in the classical Ising scenario re¬ 
quires an infinite series larger magnetic unit cells, one 
for each value of the magnetization. [5] One solution 
to the frustration problem is to avoid static order alto¬ 
gether and form a spin liquid state. |3HS] Another is to 
form ordering patterns on longer length scales via long- 
wavelength modulations or phase segregation [BHI31- Fi¬ 
nally, partial disorder, where some spins are locked into 
order and others free to exhibit Curie behavior, is seen in 
several triangular-lattice antiferromagnets. Mekata [14] 
first described the partially-disordered antiferromagnet 
(PDA) on stacked antiferromagnetic triangular lattices. 
In the PDA state, the chains of ferromagnetically-coupled 
spins form a magnetic state such that two chains are anti¬ 
aligned and one is disordered and free to flip with a low 
energy barrier. In the predicted phase diagram, the PDA 
state evolves into an ’up up down’ state with one third 
of the saturation magnetization as a function of temper¬ 
ature or magnetic held. Several stacked triangular lat¬ 


tice antiferromagnets show evidence of PDA states and 
some indication of 1 /3 plateaus including members of the 
ABX 3 and AaBB’Oe families. [51 [TUI 

Here we present the compound CoCl 2 - 2 SC(NH 2 ) 2 , de¬ 
termine its crystal structure, investigate its thermody¬ 
namic properties, and construct the magnetic phase dia¬ 
gram down to 50 mK and up to 60 T. This temperature 
and held range allows us to investigate the phase dia¬ 
gram up to saturation for magnetic helds perpendicular 
and parallel to the Ising axis, and to quantify the Ising 
anisotropy. For perpendicular magnetic helds, a trans¬ 
verse Ising scenario applies, allowing us to observe the 
held-induced quantum phase transition in a transverse 
Ising model. The structure of this material contrasts 
with the classic PDA model in that the hexagonal lay¬ 
ers in the ab plane are staggered, not stacked along the 
c-axis. Each spin in a given plane lies at the center of a 
triangle of spins in the next plane. The hexagonal lattice 
is distorted, e.g. the three exchange interactions Ji, J 2 
and J 3 in the triangle correspond to different types of 
bonds. 

In addition to the magnetic properties, we also investi¬ 
gate multiferroic behavior in this material. The coupling 
between magnetic and electric long range order is known 
as the magnetoelectric multiferroic effect in which the 
magnetic order is modified by electric field and/or the 
ferroelectricity by a magnetic field [HJ |US|. Most re¬ 
search to date in multiferroics has focused on transition- 
metal oxides. Coordination compounds are an alternate 
route to creating magnetoelectric multiferroic behavior 
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FIG. 1. a) Top view (along the c-axis) of the distorted hexag¬ 
onal planes, showing two layers, b) A CoCl2-2SC(NH2)2 
molecule c) Side view of the stacking of the planes (viewed 
along the a-axis). d) Spin level diagram of the Co^^ S = 3/2 
spin levels. 


[26l - l2^ with soft and sometimes designable lattice struc¬ 
tures. Crystallized organic molecules of thiourea [50] . 
SC(NH 2 ) 2 , and croconic acid |5T|, H 2 C 5 O 5 , are exam¬ 
ples of organic ferroelectrics. The coordination com¬ 
pound NiCl 2 - 4 SC(NH 2)2 (DTN) [28| is an example of 
a thiourea-containing compound which the polar crys¬ 
tal structure is subject to magnetostriction by magnet¬ 
ically ordered spins, creating magnetoelectric coupling. 
In the compound studied here, CoCl 2 - 2 SC(NH 2 ) 2 , two 
thiourea molecules and two Cl atoms form a tetragonal 
arrangement around each Co ion. We will show that 
magnetostrictive distortions due to either exchange in¬ 
teractions or single-ion anisotropy can modify the bulk 
electric polarization. 


II. EXPERIMENTAL 

Single crystals of CoCl 2 - 2 SC(NH 2)2 were grown at 
EWU in ethanol solution with additional isopropanol. 
Initially, C 0 CI 2 and SC(NH 2)2 were dissolved in warm 
ethanol in separate glass beakers. After mixing the so¬ 
lutions, additional isopropanol was added. After slow 
evaporation of the solvent for two week, dark blue crys¬ 
tals were obtained. Well crystallized rod-shaped crystals 
were obtained with the rod-axis oriented along the crys¬ 
tallographic c-axis. 

X-ray diffraction data were collected in the MPA-II 
group at LANL on a Bruker D 8 diffractometer, with 
APEX II charge-coupled-device (CCD) detector, and an 
American Cryoindustries Cryocool low temperature de¬ 
vice that cooled the sample to 140 K. The instrument 
was equipped with a graphite monochromatized MoK^ 
X-ray source (A = 0.71073 A), and a 0.5 mm monocap¬ 
illary. Crystals of CoCl 2 - 2 SC(NH 2)2 were mounted in a 


nylon cryoloop using Paratone-N oil. A hemisphere of 
data was collected using oj scans, with 10 -second frame 
exposures and 0.5° frame widths. Data collection and 
initial indexing and cell refinement were handled using 
APEX II software. Frame integration, including 
Lorentz-polarization corrections, and final cell param¬ 
eter calculations were carried out using SAINT- 1 - [55] 
software. The data were corrected for absorption us¬ 
ing redundant reflections and the SADABS [53] program. 
Decay of reflection intensity was not observed as moni¬ 
tored via analysis of redundant frames. The structure 
was solved using Direct methods and difference Fourier 
techniques. All hydrogen atom positions were idealized, 
and rode on the atom they were attached to. The final 
rehnement included anisotropic temperature factors on 
all non-hydrogen atoms. Structure solution, refinement, 
graphics, and creation of publication materials were per¬ 
formed using SHELXTL [55] . 


Physical property measurements were performed in the 
MPA-CMMS group at LANL, which includes the Na¬ 
tional High Magnetic Field Laboratory Pulsed Field Fa¬ 
cility. The temperature and magnetic field dependence of 
the magnetization, M(T,H), were measured in a Quan¬ 
tum Design (QD) Physical Property Measurement Sys¬ 
tem (PPMS) with a Vibrating Sample Magnetometer 
(VSM) option up to fioH = 13 T and down to T = 2 
K. Below r = 2 K, M{T,H) for H || c was measured 
with a capacitive Faraday magnetometer in an Oxford 
dilution refrigerator in a 12 T superconducting magnet. 
In addition, M{H) measurements for H |j c and H T c 
were extended up to = 60 T and down to T = 0.5 
K in resistive pulsed magnets (10 ms rise and 40 ms de¬ 
cay time) driven by a capacitor bank. The pulsed-field 
magnetization is measured via a compensated induction 
coil magnetometer [36] with in-situ sample-in sample-out 
background subtraction. Samples were immersed in ^He 
liquid or gas and their temperature was recorded by a re¬ 
sistive thermometer at zero field just prior to the pulse. 
Specihc heat, Cp{T), of sample (0.78 mg) was measured 
by the relaxation technique down to T = 50 mK in a QD 
PPMS with a dilution refrigerator option. 


The electric polarization change, AP{H), was mea¬ 
sured in pulsed magnetic fields for AP T c and for both 
H 11 c and H T c [27l|28]. Platinum contacts were sput¬ 
tered onto the samples with a cross-sectional area of 1.2 
X 1.25 mm^ for H T c and 1.75 x 2.6 mm^ for H || c. 
The induced magnetoelectric currents (analogous to py¬ 
roelectric currents) due to changes in surface charge as P 
changes with H were recorded with a Stanford Research 
570 current to voltage amplifier. Since the magnetoelec¬ 
tric current is proportional to dAP(/)/di, the measured 
signal is integrated as a function of time to obtain AP{H) 
with a high sensitivity due to the speed of the pulsed 
magnetic fields. [571HH] 
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FIG. 2. Specific heat, Cp, data of CoCl2-2SC(NH2)2, taken at 
various magnetic fields applied along H || c. Vertical arrows 
indicate two phase transition temperatures for fioH — 0.2 T. 
Inset shows the calculated total entropy, S{T), in units of 
Mn(2) by integrating C/T vs T at various magnetic fields. 


III. RESULTS 

A. Crystal Structure 

The results of single-crystal X-ray scattering and re¬ 
finement at 140 K are listed in Table HI We find the 
monoclinic structure (Cc No. 9, a = 8.199(1)A, b = 
11.542(2)1, c = 10.804(2)1, (3 = 103.587°) shown in Fig. 

The Co ions form a distorted triangular lattice in the 
ab plane that breaks spatial inversion symmetry, and al¬ 
lows for a net electric polarization. The environment im¬ 
mediately surrounding each Co ion consists of an approx¬ 
imate tetrahedron with two Cl atoms and two S atoms 
as shown in Fig. (b). The bond lengths involving the 
Co atoms are C 0 -CI 2 = 2.264(7)1, Co-Cli = 2.284(1)1, 
C 0 -S 2 = 2.309(9)1, and Co-Si = 2.324(0)1 and the an¬ 
gles around the Co atoms are CI 2 -C 0 -CI 1 = 107.88(4)° 
and S 2 -C 0 -S 1 = 96.85(4)°. 


B. Specific Heat 

The specific heat, Cp{T,H) is presented in Fig. At 
H = 0, C{T) exhibits a A-shaped anomaly at Tjv = 0.82 
K indicating a phase transition. A small upturn at the 
lowest temperatures, which is not seen in applied fields, 
may indicate a second phase transition below 50 mK. The 
anomaly centered at 0.82 K shifts to lower T in applied 
H along the c-axis and vanishes by 0.6 T. For H = 0.2 
and 0.4 T, a second anomaly is also observed as indicated 
by arrows in Fig. The entropy change occurring in the 
phase transitions was estimated by integrating CpjT as 
a function of T. The total entropy change IS.S at several 


magnetic fields is plotted in the inset of Fig. At these 
temperature, the lattice contribution is negligible and the 
total entropy change is dominated by the magnetic con¬ 
tribution. At = 0, the entropy removed by the long 
range ordering is 2.1 J/mol K, which is equivalent to 40 
% of i?ln(2). Thus the phase transitions correspond to 
only partial magnetic ordering. 


C. Magnetization vs magnetic field 

Figure [^a) shows the magnetization M{H) up to 
= 60 T in pulsed fields measured for H || c and H 
T c, calibrated to data taken in a superconducting mag¬ 
net. The data is reversible between up and down sweeps 
of the magnetic field, indicating no measurable heating 
or cooling effects during the magnetic field pulse. The 
saturation magnetization by 60 T is 3.5 /rs/f.u. for H 
II c and 2.8 /rs/f.u. for H T c, which is close to the 
expected value for S = 3/2 with minimal orbital contri¬ 
bution. The magnetic field needed to reach saturation is 
highly anisotropic, being 0.6 T for H || c and ~40 T for 
H T c. 

Given an S = 3/2 ion with a Kramer’s doublet ground 
state, there are two possible options for the ground state, 
regardless of the crystal electric field environment: an 
15”° = ±1/2) doublet ground state with an l^^ = ±3/2) 
doublet excited state, or a 15° = ±3/2) doublet ground 
state with an |5° = ±1/2) doublet excited state. Spin- 
orbit interactions can also mix the different 5° states but 
this is can be treated as a higher-order correction for a 
3d ion like Co. In our system, the high-held magnetiza¬ 
tion data is consistent only with the |5° = ±3/2) ground 
state scenario. Thus for the hard axis H ± c, saturat¬ 
ing the magnetization requires « 40 T to overcome the 
anisotropy energy and rotate the spins to the hard axis, 
while for the easy axis H || c, the magnetization satu¬ 
rates by only 0.6 T once the antiferromagnetic order is 
destroyed. On the other hand, the |5° = ±1/2) ground 
state scenario can be ruled out because it would produce 
saturation helds that are more isotropic. It cannot ac¬ 
count for our observed 60x difference between the satura¬ 
tion helds for H || c and H ± c. The largest energy scale 
in our system is the Ising anisotropy energy, and in the 
|5° = ±1/2) ground state scenario, this energy splitting 
needs to be overcome via the Zeeman effect in order to 
reach the 15^ = 3/2) for H || c but also in order to reach 
the 15a; = 3/2) state for H ± c. (|5a; = 3/2) contains a 
component of 15^ = 3/2)). Therefore we conclude that 
the Co^+ 5 = 3/2 ions have a |5° = ±3/2) ground state 
with strong Ising anisotropy and the easy axis along c. 

Moving on to the low-held behavior of M{H), this is 
shown in the inset to Fig. [^a) for pulsed helds for both 
H II and ± c, at T = 0.5 and 1.4 K, and also in Fig¬ 
ure 1^ in a superconducting magnet for H || c down 
to T = 0.1 K. The data in the different magnets are 
consistent with each other. Above = 0.8 K, M{H) 
shows Brillouin-like behavior. Below T/y, M{H) forms a 
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FIG. 3. a) High-field M(H) for H || c and H _L c at T = 0.5 and 1.4 K, taken in pulsed magnetic fields up to 60 T. Vertical 
arrow near H = 30 T indicates a slope change in dM{H)/dH. Inset shows an expanded scale at low fields, b) Magnetization 
isotherms, M{H), of CoCl2-2SC(NH2)2 for H || c at selected temperatures. Inset shows M{H) up to 12 T taken at T = 0.2 K. 
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FIG. 4. Inverse magnetic susceptibility, H/M{T), of G 0 GI 2 - 
2 SG(NH 2)2 as a function of temperature for two orientations 
of the magnetic field; H || c and H _L c. Solid lines represent a 
Gurie-Weiss fit. Inset shows magnetic susceptibility below T 
= 1.5 K at selected magnetic fields applied along H || c. Phase 
transition temperatures, determined by d[{M/H) ■ T]/dT, are 
indicated by arrows. 


plateau at l/3Msat for H\\c, and saturates for fj^oH > 0.6 
T with a saturated magnetization of 3.2 /rs/f.u. The 
onset in magnetic field of the 1/3 step, and of the sat¬ 
uration correspond to peaks in the heat capacity shown 
previously. 


D. Magnetic susceptibility 


The inverse magnetic susceptibility, H/M{T), is plot¬ 
ted in Fig. 1^ A fit by a Curie-Weiss law to the data 
above T = 100 K results in a Curie temperature 9p = 11 
K and an effective moment fies = 4.4 /is for H || c and 
= —15 K and = 4.1 for H T c, respectively. 
The moments are in agreement with the expected Co^"*" 
ion value with spin S = 3/2. The large anisotropy in the 
Curie-Weiss temperature, with a sign change from H || 
c to H T c, indicates that the origin of the Curie-Weiss 
behavior is primarily the single-ion anisotropy while the 
exchange interactions play a smaller role |37j . This is con¬ 
sistent with the low values of T^r = 0.8 K and He = 0.6 
T relative to 9p. 


The low-temperature M{T)/H curves at selected mag¬ 
netic fields are plotted in the inset of Fig. The 
M{T)/H curve at H = 0.05 T reveals a kink at Tn = 
0.83 K, below which M{T)/H diverges with decreasing 
temperature. As the magnetic field increases, Tjv shifts 
to lower temperatures and vanishes for fJ^gH > 0.6 T. 
A second anomaly in M{T)/H curves is also seen for 
0.1 < < 0.4 T (see the representative = 0.2 

T curve in Fig. |^. This anomaly corresponds to the 
second phase transitions previously shown in the heat 
capacity that is the onset of the ”1/3” state with 1/3 of 
the saturation magnetization. Fig. shows that inside 
the 1/3 state, the magnetic susceptibility approaches a 
constant value as T —)■ 0, rather than diverging as it does 
at H = d. 
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FIG. 5. (a) Electric polarization change with time, dAP/dt, plotted as a function of magnetic field ^oH at various temperatures, 
where H || c and H _L AP. Inset shows the dAP/dt at T = 0.48 K up to P = 60 T. (b) AP{H) determined by integrating the 
dAP/dt as a function of time. 




FIG. 6. (a) Electric polarization change, dAP/dt vs magnetic field noH, at T = 0.5 K for H _L c. Inset shows an expanded 
plot up to fj,oH = 10 T. The vertical arrow indicates a phase transition field and horizontal arrows indicate up-sweep and 
down-sweep, (b) AP{H) determined at various temperatures by integrating the dAP/dt as a function of time. Inset shows 
the dAP/dt at selected temperatures. 


E. Electric Polarization 

The electric polarization change, AP{H), for P T c as 
a function of magnetic held was measured in pulsed mag¬ 
netic Helds up to 60 T. Consistent with the intrinsically 
polar crystal structure, no dependence on applied voltage 
up to 200 V is seen before or during measurements. The 
measurement of magnetoelectric current for H || c and H 
T AP was performed in a slow capacitor-driven pulsed 
Held magnet (40 ms rising and 250 ms falling time) up to 
5 T, and in a faster capacitor-driven pulsed magnet (10 
ms rising and 40 ms falling time) up to 60 T. The mea¬ 
sured magnetoelectric current, d(AP) /dt, induced by the 
polarization change of the sample and the integrated sig¬ 
nal, AP{H), are plotted in Fig. [^(a) and (b) for H || c 
and in Fig. |^(a) and (b) for H T c, respectively, where 
AP T c for both orientations of the magnetic field. The 
difierence in magnitude of the raw d(AP)/dt data be¬ 


tween Fig. (a) and the inset is due to the difierent 
magnetic field sweep rates. This sweep-rate dependence 
is absent for the integrated AP{H) curves. 

AP{H) in Fig. [^shows features consistent with M{H) 
for both magnetic field directions. For H || c, a 1/3 
plateau is seen in P{H) similar to M{H), at T = 0.45 K. 
For T > Tn, AP{H) evolves smoothly and monotoni- 
cally without any noticeable anomalies and the magni¬ 
tude is immediately suppressed above Tjy. For H T c, a 
kink in d(AP)/dt appears at fioH = 1.5 T correspond¬ 
ing to the suppression of long-range order (see inset) and 
a drastic change occurs near hqH = 30 T. The inte¬ 
grated P{H) shows a large peak at 30 T, which is just 
below the field where the magnetization begins to sat¬ 
urate. The same behavior is seen for both rising and 
falling field sweeps, except for a small hysteresis at the 
phase transition near (J-qH = 1.5 T. The observed ampli¬ 
tude of AP{H) at the phase transition is similar for both 
























6 


H II c and H _L c, whereas AP{H) around fXoH = 30 T 
along H _L c is ~ 100 times bigger. 


IV. DISCUSSION 

The magnetic ordering in CoCl 2 - 2 SC(NH 2)2 for H |j c 
is bounded by T = 0.83 K and H = 0.6 T. The phase 
diagram obtained by the present measurements is shown 
in Fig. (a), where the phase boundaries were deter¬ 
mined from the peak position in d[{M{T)/H) ■ T]/dT, 
dM{H)/dH, d[AP{H)]/dH, and the peak position of 
Cp{T). The phase diagram and magnetic behavior is 
consistent with a PDA model. m In the zero-field par¬ 
tially disordered state, the magnetization diverges at low 
temperatures, indicating some free spins that exhibit a 
Curie-behavior. In this phase, only ^ 40% of the ex¬ 
pected entropy for fully ordering the spins is observed 
from the specific heat measurements. Upon applying a 
field, a phase transition to the 1/3 magnetization phase 
is observed in the heat capacity, magnetization, and elec¬ 
tric polarization measurements. In this phase with 1/3 
of the saturation magnetization, the Curie-like magnetic 
susceptibility disappears and the magnetization at low 
temperatures approaches a constant value -all spins lock 
into the 1/3 state. The 1/3 magnetization holds over 
an extended region of the T — iJ phase diagram with no 
observable change in magnetization with temperature or 
field. This behavior is an especially pronounced exam¬ 
ple of the behavior expected from Ising spins in a PDA 
model. The magnetic field required to order the disor¬ 
dered spins and induce the 1/3 state extrapolates close 
to i/ = 0 at T = 0 , indicating minimal energy barriers 
for flipping the free spins. Given Ising spins on a lattice 
with 3-fold symmetry (approximately satisfied for this 
compound), an ’up up down’ configuration within trian¬ 
gles is the most likely ordering. With further increasing 
field the spins undergo a second held-induced phase tran¬ 
sition out of this locked 1/3 state and then evolve towards 
saturation at 3 /is. 

An important difference between this material and 
other PDA materials is a lack of 1-D chains in the crys¬ 
tal structure. The classic PDA model [T3] is postulated 
for strongly-correlated ferromagnetic chains of spins that 
are coordinated in a triangular motif in the perpendicu¬ 
lar direction. In CoCl2-2SC(NH2)2, the crystal structure 
does not evidence such chains - the hexagonal layers are 
staggered, not stacked. A spin in one plane is bonded 
to three spins in the next plane, forming a tetrahedron. 
Thus it is a geometrically different example of partial 
antiferromagnetic ordering due to triangular frustration. 

[siiiniiiiiinHis] 

We can quantify the degree of Ising-ness of the Co 
spins from the anisotropy between the saturation mag¬ 
netic held for the easy c-axis (0.6 T) and the hard axis 
perpendicular to c (40 T). The scenario consistent with 
the data is for the Co^’*' S = 3/2 spins to be split into 
an |S'^ = ±3/2) Kramer’s protected ground state, with 


an excited = ± 1 / 2 ) state separated by an energy gap 
D (Fig. (d)). For H || c, saturation is reached for 
low helds ^s soon as effects of magnetic exchange inter¬ 
actions J are overcome by 0.6 T) while saturating the 
magnetization for H ± c requires mixing components of 
| 5'2 = ±3/2) and = ±1/2) to achieve the \Sx = 3/2) 
state. Thus the gap D must be closed to reach H ± c 
saturation near 40 T, and so we estimate D ~ 60 K along 
the c-axis, assuming g = 2 and S = 3/2. 

A minimal spin Hamiltonian that captures magnetic 
exchange and single-ion anisotropy can be written as 
H = Er,v ■ Sr+e„ ±E.[^(^/)'-^ 7 MBS^H], where 

Jv are the magnetic exchange constants, D is the single¬ 
ion anisotropy due to a uniaxial crystal electric held, and 
By = (ax, by, xz) are the relative vectors between nearest- 
neighbor Co-ions connected by magnetic exchange inter¬ 
actions. From the crystal structure we know that the 
uniaxial anisotropy axis lies along the c-axis for both Co 
sites in this material, and the magnetic properties show 
this to be the easy axis. 

A comparison of the saturation helds of 0.6 T along 
the easy axis (needed to overcome antiferromagnetic ex¬ 
change interactions) and 40 T along the hard axis (needed 
to overcome single-ion anisotropy) indicates that this ma¬ 
terial is in the limit of large D/J with an Ising energy 
scale approximately 10 times larger than the magnetic 
exchange. A consequence of the large D/J limit is that 
the Curie-Weiss behavior of the magnetic susceptibility is 
dominated by single-ion anisotropy and not by magnetic 
exchange. m Consistent with this, the Curie-Weiss tem¬ 
perature changes sign when the magnetic held is rotated 
from the easy to hard axis (-15 K vs -1-11 K). 

For the magnetic held H ± c, a transverse Ising sce¬ 
nario applies. The magnetization shows three diherent 
slopes, between 0<H<2T, 2<H<40T, and 
40 < i? < 60 T. For lower helds, we can expect the 
magnetization to be inhuenced by spin hip excitations 
(superposition states of jS"* = 3/2) and \S^ = —3/2)). A 
phase transition near 2 T at 0.5 K in M(H) is observed 
in Fig. which would be a held-induced quantum phase 
transition of the transverse Ising model. [38114(1] Further 
investigation of transverse Ising behavior in this system 
is a future project requiring a detailed understanding of 
the exchange couplings. The hnal saturation near 60 T 
is given by the energy scale needed to overcome D, the 
gap between the ground state = ±3/2) and excited 
1 = ± 1 / 2 ) states. 

Coupling between magnetism and electric polarization 
is shown in Fig. (b). CoCl 2 - 2 SC(NH 2)2 has a polar 
crystal structure that allows for a net electric polariza¬ 
tion. Thus it is a type of ferroelectric with an order¬ 
ing temperature above the melting point of the crystal. 
Any magnetostriction-induced changes in the lattice pa¬ 
rameters can therefore modify the electric polarization. 
In general, the quantities D and J depend on the ar¬ 
rangement of atoms within the crystal structure. Thus, 
the system will modify these through magnetostriction 
to minimize the combination of magnetic and crystalline 
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FIG. 7. The H -T phase diagram of CoCl2-2SC(NH2)2 for 
H II c, obtained from magnetization, specific heat and elec¬ 
tric polarization measurements. Phases are indicated as para¬ 
magnetic (PM), partially-disordered antiferromagnet (PDA), 
l/3Maat, and Msat for the phase approaching saturation, (b) 
The electric polarization change, AP{H) left axis, at T = 0.5 
K and magnetization, M{H) right axis, at T = 0.1 K, which 
are the respective base temperatures of measurements. 


energy. In the data, we observe that in the region of mag¬ 
netic ordering, AP{H) exhibits similar magnetic field 
evolution as the magnetization including a 1/3 plateau 
and paramagnetic behavior just above T/y. The magni¬ 
tude of AP is 3 /iC/m^ for H || c due to magnetic ex¬ 
change striction. Much larger values of AP{H) are seen 
for H _L c at higher magnetic fields up to 60 T. We can 
attribute these to magnetostriction driven by the single¬ 
ion D term in the Hamiltonian. The value of AP{H) is 
300 which is within a factor of ten of the largest 

magnetic field-induced electric polarization changes ob¬ 
served in multiferroic materials mma. The thiourea- 
containing compound NiCl2-4SC(NH2)2 also shows mag¬ 
netoelectric coupling due to both D and J terms, albeit 
with smaller magnitude. A detailed experimental and 
theoretical analysis of these effects are presented in Ref. 



V. SUMMARY 

Single crystals of a new coordination compound C 0 CI 2 - 
2 SC(NH2)2 have been synthesized. The compound has a 
distorted hexagonal structure and shows magnetic order¬ 
ing consistent with a partially-disordered antiferromag¬ 
netic state at H = 0 where two spins order and one re¬ 
mains disordered in a triangular motif. In applied mag¬ 
netic fields, a locked 1/3 state occurs, where the magneti¬ 
zation remains constant at 1 /3 of the saturation magne¬ 
tization over a remarkably broad region in temperature 
and magnetic fields. The geometry of this compound dif¬ 
fers from that of usual partially disordered antiferromag- 
nets, since it shows staggered hexagonal planes, rather 
than the usual c-axis spin chains in a triangular configu¬ 
ration. The 00^+ S = 3/2 spins form a doublet ground 
state with 15^ = 3/2) and an easy-axis anisotropy axis 
along the crystallographic c-axis. The well-separated en¬ 
ergy scales for magnetic ordering and anisotropy make 
this a clean Ising system. In addition, this material 
shows an electric polarization that is strongly coupled to 
magnetization with a magnetic field-induced polarization 
change up to 300 /rC/m^. The largest electric polariza¬ 
tion changes in this material are produced by single-ion 
anisotropy effect creating magnetically-driven distortions 
of the Co environment. 
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TABLE I. Crystallographic data 


Formula 

CoCl2-2SC(NH2)2 

Formula Weight 

282.07 

T(K) of measurement 

140 

Crystal system 

monoclinic 

space group 

Cc 

a (A) 

8.1991(15) 

6(A) 

11.542(2) 

c (A) 

10.804(2) 

a(“) 

90 

pn 

103.587(3) 

7n 

90 

Volume (A®) 

993.9(3) 

Z 

4 

Calculated p(mg/m^) 

1.885 

Crystal Size(mm®) 

0.20x0.14x0.14 

Transmission Tmi-a, Tmax 

0.621, 0.710 

Data Collection 6min, Omax (°) 

3.11, 28.44 

R (reflections) 

0.0309 (2274) 

wR2 (reflections) 

0.074 (2308) 

GOF of 

1.166 

Ap (max., min.) (eA“®) 

1.207, -0.891 



